Nonlinear electrohydrodynamics of a viscous droplet 
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A classic result due to G.I.Taylor is that a drop placed in a uniform electric field becomes a 
prolate or oblate spheroid, which is axisymmetrically aligned with the applied field. We report an 
instability and symmetry-breaking transition to obliquely oriented, steady and unsteady shapes in 
strong fields. Our experiments reveal novel droplet behaviors such as tumbling, shape oscillations, 
and chaotic dynamics even under creeping flow conditions. A theoretical model, which includes 
anisotropy in the polarization relaxation due to drop asphericity and charge convection due to drop 
fluidity, elucidates the interplay of interfacial flow and charging as the source of the rich nonlinear 
dynamics. 
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Nonlinear phenomena such as instabilities and turbu- 
lence naturally occur in fluid dynamics because of the 
nonlinearity of the Navier-Stokes due to the inertial term. 
In the absence of inertia, the Stokes equations governing 
the fluid flow are linear and evolving boundary conditions 
are the only source of nonlinearity pQ. For example, a sin- 
gle particle exhibits complex dynamics if it is deformable: 
a capsule or a red blood cell in shear flow [MS]) a drop 
in oscillatory strain-dominated linear flows [7] , or a drop 
sedimenting in an electric field [5]. 

If the particle is rigid, chaotic motions are observed 
with ellipsoids in shear flow [9] or uniform electric fields 
[T(3| ITT] . In the latter case, the nonlinear dynamics arises 
from anisotropy in the polarization relaxation due to as- 
phericity of the particle shape. A question arises - what 
if the particle is not rigid, but fluid and its shape is not 
fixed? Would particle fluidity and deformability give rise 
to new, richer behaviors? In this Letter, we explore these 
questions on the example of a viscous drop subjected to a 
uniform direct current (DC) electric field. Upon increase 
in the field strength, this system was found to undergo a 
symmetry-breaking transition from axisymmetric to lin- 
ear flow illustrated in Figure [l] B and C [12], which 
hints upon the possibility of more complex physics. We 
first discuss nonlinear drop electrohydrodynamics theo- 
retically, and then we describe its experimental realiza- 
tion. We find that depending on the fluids viscosity, the 
droplet may exhibit a range of intriguing dynamics such 
as tumbling, sustained shape oscillations, and even chaos. 

Theoretical model for drop polarization and shape evo- 
lution. When placed in an electric field, a particle po- 
larizes because free charges carried by conduction accu- 
mulate at boundaries that separate media with different 
electric properties. This is illustrated in Figure [2] on the 
example of a sphere in a uniform electric field. The mag- 
nitude and orientation of the induced dipole depend on 
the mismatch of electric properties between the particle 
( "in" ) and the suspending fluid ( "ex" ) 
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FIG. 1: Exposure to a uniform direct current (DC) electric field 
with increasing strength excites a variety of drop responses |12| 
(see videos \13\). A. A castor oil drop suspended in silicon oil is 
spherical in the absence of electric field. B. Weak fields induce 
axisymmetric oblate deformation. C. In strong fields, the drop is 
tilted with respect to the applied field direction and the flow has a 
rotational component. The sketches illustrate the drop shape and 
flow streamlines. 

where a and e denote conductivity and dielectric con- 
stant, respectively. The product of R and S compares 
the charge relaxation times of the media [TU [TS] 

RS = , where t cM = — r c , ox = — . (2) 
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If RS < 1 (r c i n > r c cx ), the conduction response of 
the exterior fluid is faster than the particle material. As 
a result, the induced dipole is oriented opposite to the 
applied electric field direction. This configuration is un- 
favorable and becomes unstable above a critical strength 
of the electric field [TU1 HE1 HZ] ■ A perturbation in the 
dipole alignment gives rise to a torque, which drives phys- 
ical rotation of the sphere. The induced surface-charge 
distribution rotates with the particle, but at the same 
time the exterior fluid recharges the interface. The bal- 
ance between charge convection by rotation and supply 
by conduction from the bulk results in an oblique dipole 
orientation shown in Figure [2j(b). The rotation rate ui 
is determined from the balance of electric and viscous 
torques acting on the particle, P x E = A • a>, where A 
is the friction matrix. 

The spontaneous spinning of a rigid sphere in a uni- 



FIG. 2: Charge distribution and induced dipole P for a sphere 
with RS < 1. Above a critical field strength Eq > Eq , where 
Eq = |E| and Eq is given by Eq. HI constant rotation around an 
axis perpendicular to the electric field is induced by the misaligned 
dipole moment of the particle (right). The rotation can be either 
clock- or counter-clockwise. 



FIG. 3: Sketch of the model for drop electrorotation. Drop as- 
phcricity is characterized by the aspect ratio f3=a«/a±. The rigid 
body rotation is related to the tilt angle of the deformed shape, 
(9, by u> = d6/dt. In the co-rotating frame the interface velocity is 



form DC electric field has been known for over a cen- 
tury, first attributed to the work of Quincke in 1896. In 
this case, the friction matrix is diagonal, and a straight- 
forward calculation |16H18j . assuming instantaneous po- 
larization, yields three possible solutions: no rotation, 
u> = 0, and 
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where the ± sign reflects the two possible directions of 
rotation and 
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Tmw, the Maxwell- Wagner polarization time, is the char- 
acteristic time-scale for polarization relaxation. The 
dipole "tilt" is steady; the angle between the dipole and 
the electric field is (f> = arctan [(T mto aj) -1 ] . Eq. [3] shows 
that rotation is possible only if the electric field exceeds a 
critical value given by Eq. Hence, if Eq < Eq, the sphere 
and the suspending fluid are motionless. If Eq > Eq, the 
sphere rotates and drags the fluid in motion; the resulting 
flow is purely rotational. 

If polarization relaxation, i.e., non- instantaneous 
charging of the interface described by <9 t P = uj x P — 
T mw (P — Peg), is included in the analysis, the polariza- 
tion evolution equation and the torque balance map onto 
the Lorenz chaos equations QUI [H] • A second bifurcation 
occurs in stronger fields and the sphere exhibits chaotic 
reversal of the rotation direction. 

Unlike solid particles, drops are fluid and have a free 
boundary. The electric stress deforms the drop and, as 
a result, the friction matrix and polarization relaxation 
become anisotropic and dependent on the drop orienta- 
tion relative to the applied electric field. Moreover, the 
interface does not move as a rigid body and its veloc- 
ity differs from that the Quincke rotation. The modified 
torque balance (assuming for simplicity that drop shape 
remains an axisymmetric oblate spheroid with fixed as- 
pect ratio once rotation is initiated) is 
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where a± is the friction coefficient, x°° is the high- 
frequency susceptibility, and ||, _L denote components 
parallel and perpendicular to the axis of symmetry, see 
Figure [3] Drop fluidity is accounted by interfacial ve- 
locity with frequency v. The latter is determined from 
the balance of electrical energy input and viscous dissipa- 
tion, in a similar fashion to the analysis of drop rotation 
in a magnetic field [20 , or red blood cell tank-treading 
in shear flow |21) . The polarization relaxation equations 
in a coordinate system rotating with u> are 
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where ry^ and x» ± ar e directional Maxwell- Wagner re- 
laxation timescales and low-frequency susceptibility re- 
spectively. If v = 0, Eq. [5] and Eq. [6] reduce to the 
equations of motion for a rigid ellipsoid [IT] , and predict 
three types of behavior: alignment of the symmetry axis 
with the electric field, oscillations around the field direc- 
tion ("swinging"), and continuous flipping ("tumbling"). 
The additional torque associated with the interface flu- 
idity (characterized by v) modifies these behaviors to: 
steady axisymmetric orientation (Taylor regime), steady 
tilted orientation, swinging around a non-zero tilt angle 
with respect to the electric field, and tumbling. More- 
over, chaotic switching between the swinging and tum- 
bling states is also found. The phase diagram resulting 
from the numerical solution of Eq. [5] and Eq. [6] is shown 
in Figure [4] In all cases except for the chaotic behavior, 
the motion of the dipole mirrors the motion of the el- 
lipsoid. The model predicts that chaos is observed with 
high viscosity drops (drop viscosity at least four times 
greater than the suspending medium). In this case, the 
chaotic motion stems from the lack of synchronization 
between ellipsoid and dipole orientations. 
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FIG. 4: Phase diagram for a fluid ellipsoid with viscosity ratio A = 
Min/^ex = 14. The lines represent the boundaries between various 
behaviors computed from the model. The symbols correspond to 
experimental data: □ - Taylor (steady axisymmetrically oriented 
ellipsoid), o - steady tilted ellipsoid, A - chaotic tumbling. The 
lines crossing the symbols reflect the variations in the aspect ratio 
of the drop during the rotation. 

Of course in reality drop shape does not remain an el- 
lipsoid with fixed aspect ratio, and the experiments show 
pronounced shape oscillations in the case of low- viscosity 
drops (see Figure [5] a). This is due to the fact that the 
rotation period is comparable with the surface tension 
relaxation time a/i cx /7 (7 being the interracial tension, 
//ex the suspending fluid medium, and a the initial drop 
radius). However, in this case the main axis is found to 
oscillate around a defined tilted angle and hence this be- 
havior qualitatively corresponds to the swinging mode. 
High- viscosity drops remain nearly undeformed through- 
out one period of rotation and tumble like rigid ellipsoids 
(see Figure [5jb) 
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FIG. 5: Examples of oscillatory and tumbling drop behavior, 
a) A=l, S =9.9kV/cm, a=1.8mm. b) A=14, £ =9.7kV/cm, 
a=3.0mm 

Our "fluid ellipsoid" model for drop electrorotation 
identifies drop viscosity as a key control parameter for 
drop behavior. Next we experimentally test the model 
and show that it qualitatively captures the variety of drop 
responses. 

Experiment: Weakly conducting fluids are used for the 
drop (silicone oil) and continuous phase (castor oil) . This 
experimental system is characterized by conductivity ra- 
tio R = 0.03 and permittivity ratio S — 1.8. To in- 



vestigate the effect of drop viscosity, the viscosity ratio 
A = /Uin/j^ex is varied in the range 1 to 14. The electric 
field is increased in small steps of about O.IEq and the 
system is allowed to equilibrate at each step to avoid spu- 
rious transients. For example, it is possible to excite drop 
tumbling upon a sudden increase of the field strength but 
the drop eventually settles into steady tilted state. More 
details about the experiment can be found in [12]. 

Typical drop behaviors are illustrated in Figure [6] The 
unsteady dynamics of high viscosity drops (A = 14) is 
chaotic tumbling, while low viscosity drops (A = 1) un- 
dergo steady shape oscillations. Intermediate viscosity 
drops (A = 4) exhibit both behaviors, namely a cycle 
consisting of shape oscillations with increasing amplitude 
followed by several tumbles. These behaviors are better 
seen in the insets of Figure [6j where the angle between 
the drop major axis and the applied field direction is 
plotted. 

For high viscosity drops A = 14, unsteady behavior is 
observed at aspect ratio above (3 > 1.2 in both the model 
and experiment. Experimentally, we observe a perturba- 
tion to the oblate drop shape to slowly grow until reach- 
ing a chaotic tumbling state, in which the direction of ro- 
tation reverses chaotically. The transitions between the 
Taylor, tilted and tumbling states of a drop are given by 
lines in Figure |6]a. The largest stretching occurs while 
the short axis of the oblate drop is aligned with the field 
(see Figure[5jb at 0.7s and 2.37s). At this point, the drop 
can either break up or continue to rotate. 

Low viscosity drops ( A = 1) exhibit steady oscillations 
between a prolate and a tilted oblate shapes, see Figure 
[5] a. These shape oscillations can be explained by con- 
sidering the dynamics of the induced dipole and shape. 
The lower fluid viscosity allows the drop to be easily de- 
formed by the electric field acting on the induced surface 
charge. As the interface charges, the drop is squeezed 
into an oblate shape. The accumulated surface charge ro- 
tates due to the electric torque until the induced dipole 
is nearly aligned with the field. In this dipole configu- 
ration, the drop is pulled into a prolate shape until the 
surface charge relaxes. The cycle repeats by the interface 
recharging and deforming into an oblate shape. At fields 
strengths well above the transition to unsteady behavior, 
the amplitude and frequency of oscillation increase until 
the drop stretches into a highly deformed prolate " dumb- 
bell" shape, which then either undergo breakup or col- 
lapses back into itself, disrupting its rotational velocity. 
After collapsing, the oscillations will begin to build again 
in the same direction as before. 

The different unsteady behaviors are attributed to the 
interplay of shape and dipole dynamics. The more vis- 
cous the drop fluid, the higher resistance to fluid motion. 
Accordingly, high viscosity drop tend to behave as rigid 
particles and tumble in the electric field. They deform 
very little during the dipole rotation. In contrast, low 
viscosity drops undergo large deformations seen as shape 
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FIG. 6: Phase diagrams for viscosity ratios A =1, 4, 14. The elec- 
tric field is slowly increased and steady-state behavior is observed. 
Taylor (Ta) indicates axisymmetric flow and oblate deformation. 
Tilted (Ti) indicates a tilted drop orientation with rotational flow. 
Unsteady (Un) indicates time dependent drop shape and orienta- 
tion. Breakup (Br) indicates regions where drop breakup is ob- 
served. The inset shows the time dependent behavior for one drop 
indicated on the chart with o, radial units of aspect ratio j3 and 
angle units of drop orientation 9. 



oscillations. The limited deformation of high viscosity- 
drops allows for the fixed shape model to accurately pre- 
dict the transition to unsteady behavior. 

In conclusion, in this Letter we report novel nonlin- 
ear dynamics of a droplet in uniform DC electric fields. 
Droplet tumbling, shape oscillations, and chaotic rota- 
tions occur under creeping flow conditions, where non- 
linear phenomena are rare. We have constructed a phase 
digram of drop behaviors as a function of viscosity ratio 
and field strength. The experimental data qualitatively 
agrees with a theory which models the drop as a fluid el- 
lipsoid. The favorable comparison is surprising given the 
many simplifications in the theory and encouraging for 
future efforts to build a more refined and accurate model. 
The nonlinear drop physics could inspire new approaches 
in electromanipulation for small-scale fluid and particle 
motion in microfluidic technologies. 
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